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ABSTRACT 

We calculate the expected amplitude of the dipole and higher spherical harmonics 
in the angular distribution of radio galaxies. The median redshift of radio sources 
in existing catalogues is z ~ 1, which allows us to study large scale structure on 
scales between those accessible to present optical and infrared surveys, and that of 
the Cosmic Microwave Background (CMB). The dipole is due to 2 effects which turn 
out to be of comparable magnitude: (i) our motion with respect to the CMB, and 
(ii) large scale structure, parameterised here by a family of Cold Dark Matter power- 
spectra. We make specific predictions for the Green Bank 1987 (87GB) and Parkes- 
MIT-NRAO (PMN) catalogues, which in our combined catalogue include ~ 40, 000 
sources brighter than 50 mJy at 4.85 GHz, over about 70% of the sky. For these 
relatively sparse catalogues both the motion and large scale structure dipole effects 
are expected to be smaller than the Poisson shot-noise. However, we detect dipole and 
higher harmonics in the combined 87GB — PMN raw catalogue which are far larger 
than expected. We attribute this to a 2% flux mismatch between the two catalogues. 
Ad-hoc corrections to match the catalogues may suggest a marginal detection of a 
dipole. To detect a dipole and higher harmonics unambiguously, a catalogue with full 
sky coverage and ~ 10 6 sources is required. We also investigate the existence and 
extent of the Supergalactic Plane in the above catalogues. In a strip of ±10° of the 
standard Supergalactic equator, we find a 3cr detection in PMN raw , but only la in 
87GB raw . We briefly discuss the implications of on-going surveys such as FIRST and 
NVSS and follow up redshift surveys. 

Key words: galaxies: large scale structure - radio galaxies - dipole - spherical 
harmonics - etc. 



1 INTRODUCTION 

Our local universe is overwhelmingly clumpy. As we look on 
increasingly larger scales, we continue to see clustering. Ob- 
servations of discrete objects and their clustering properties 
on ~ 100/i _1 Mpc scales start to show evidence for homo- 
geneity. The largest scales we observe (~ 1000 h~ Mpc) are 
by looking at the Cosmic Microwave Background (CMB). 
In the CMB (apart from a dipole component attributed to 
our motion) we see a striking level of isotropy, with fluctua- 
tions at the level of 10 -5 on 10° scales (Smoot et al. 1992). 
Somewhere between the scale accessible to us by the CMB 
and that of galaxies and quasars, the universe changed from 
smooth to lumpy, forming structures from some small initial 
fluctuations. 

It has been suggested (Kaiser 1984) that galaxies form 
preferentially in high density peaks of the underlying mass 
distribution. If this is true, then the statistics of galaxy dis- 



tributions provide us with information, albeit biased, about 
the underlying dark matter distribution. 

Optical and infrared galaxy surveys have been used ex- 
tensively to study clustering out to ~ 200 ft -1 Mpc, or a 
redshift of z ~ .07 (e.g. Strauss & Willick 1995 for review). 
The high luminosity of radio galaxies and quasars, makes 
them detectable out to large, cosmological distances (z ~ 4) 
and consequently, useful in studying clustering out to high 
redshift. Unfortunately, most large surveys to date provide 
only fluxes and angular coordinates on the sky, not redshifts, 
for individual objects. This means that at present, only the 
projected distributions of radio galaxies can be studied. 

It has long been debated whether radio galaxies are 
clustered or isotropic on the largest scales. The much quoted 
study by Webster (1976), which looked at 8000 radio sources 
at 408 MHz, found < 3% variability in the number of sources 
found in a randomly placed 1 Gpc cube. This led to the gen- 
erally accepted view that radio sources were isotropically 
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distributed. Even if this were not true, the large range in 
intrinsic luminosities of radio sources would effectively wash 
out structures when the distribution was projected onto the 
sky and information about the objects' radial distribution 
was lost. More recently, Shaver & Pierre (1989) reported a 
detection of slight clustering of bright, nearby radio sources 
to the Super Galactic Plane (SGP) and Benn & Wall (1995) 
discussed other measures of anisotropy in radio surveys. In 
Section 7 we investigate the presence of the SGP in the 
87GB raw and PMN raw and report on marginal detection. 
Clustering in the above catalogues was also studied by cor- 
relation function analysis (Kooiman, Burns & Klypin 1995; 
Sicotte 1995; Loan, Wall & Lahav 1997). These studies in- 
dicated that radio galaxies are actually more strongly clus- 
tered than local optical galaxies. This conclusion is also con- 
firmed by correlation analysis of the FIRST survey (Cress 
et al. 1997). Motivated by these recent indicators of clus- 
tering in radio surveys and by the upcoming deeper surveys 
such as FIRST and NVSS (see Condon 1997 for review), we 
attempt here to study large scale structure in radio surveys 
using Spherical Harmonic Analysis (SHA). 

At this point, we introduce the labels that we use to 
refer to various forms of the 87GB and PMN catalogues. 
From here on, we will distinguish between the raw, uncor- 
rected, combined catalogue (87GB — PMN raw ), and our 
Af matched (corrected) version of the combined catalogue 
(87GB — PMN ma tch)- When talking about predictions, wc 
use the label 87GB-PMN (with no subscripts) to refer to an 
ideal 4.85 GHz survey. 

In Section 2, we present the 87GB raw and PMN raw 
catalogues and discuss their mismatch. In Section 3, we show 
the formalism of predicted rms harmonics for the models. In 
Section 4, we look at the radio dipole and how it compares 
to the predicted dipoles due to shot noise, our motion with 
respect to the CMB, and large scale structure. In Section 5, 
we consider the observed dipole in the combined catalogue, 
and in Section 6, apply the SHA to the catalogue. In Section 
7, we discuss the Supergalactic Plane and in Section 8, wc 
summarize the results and discuss future surveys. 



2 THE RADIO CATALOGUES 

2.1 The Green Bank survey (87GB) 

The 87GB survey was carried out by Condon, Broderick & 
Seielstad (1989), over a period of three years (1986-88), us- 
ing the NRAO 91m transit radio telescope (in Green Bank, 
West Virginia, USA) in conjunction with the NRAO 4.85 
GHz, 7-beam receiver that was built for the survey. A large 
part of the sky visible from the Northern Hemisphere was 
observed ( 0° < <5 < +75°, h < a < 24 h ). Errors due to in- 
strument noise and source confusion for the survey were esti- 
mated to be a ~ 5mJy and a detection threshold of S > 5a 
was adopted in compiling the source lists. A catalogue of 
54,579 radio sources above a flux limit of ~ 25mJy was sub- 
sequently produced (Gregory & Condon, 1991) from obser- 
vations taken in 1987. Excluding a few small regions thought 
to be contaminated by sidelobes of strong local sources or by 
solar interference, data from all the surveyed regions were 
included in the source catalogue, which covers ~ 6 steradi- 
ans on the sky. 



2.2 The Parkes-MIT-NRAO survey (PMN) 

The PMN survey was carried out by Griffith & Wright 
(1993) over the course of 1990, using the Parkes 64m tele- 
scope (near Parkes, New South Wales, Australia) in con- 
junction with the same 7-beam receiver as used in 87GB 
survey. The area to be surveyed, which covered —87.5° < 
S < +10.0°, for all right ascensions (0 h < a < 24 h ) in the 
Southern Celestial Hemisphere, was divided and observed in 
4 declination strips (see Table 1, Griffith & Wright 1993). 
The source detection reliability criterion was 90%, corre- 
sponding to about S > 4.4a. The definition adopted of a 
real source, was that its observed flux was close to the flux 
limit - not necessarily above the flux limit. Source catalogues 
for three of the four declination strips have been complied 
(Tropical - by Griffith et al. 1994; Southern - by Wright et al. 
1994 and independently by Gregory et al. 1994; Equatorial 
- by Griffith et al. 1995), together covering ~ 6.4 steradians 
of sky. More recently, the source catalogue for the fourth 
declination strip (Zenith) has become available (Wright et 
al.1996), but because of its lower flux limit of about 72mJy, 
it is not included in this study. 

2.3 How well are the resulting catalogues 
matched? 

We intend to use the 87GB raw and PMN raw catalogues 
together to provide extensive sky coverage, and need some 
idea of what inconsistencies might arise from the fact that 
different telescopes were used to collect the data and differ- 
ent reduction algorithms and selection criteria were used in 
compiling the catalogues. Fortunately, there are two regions 
of the sky for which there is overlap information. 

2.3.1 Region 1 

The declination band 0° < 8 < 10° was observed with both 
the NRAO and Parkes telescopes and the data collected re- 
duced by 87GB catalogue method (Gregory & Condon 1991; 
Gregory et. al 1994) as well as the PMN catalogue method 
(Griffith et al 1994; Griffith et al. 1995). Griffith et al.(1995) 
examine the agreement between resulting catalogues. They 
find some discrepancies in positions and a 2% disagreement 
of flux scales. A number of explanations for these inconsis- 
tencies are offered: 

• Two different telescopes were used to collect the data. 

• The observations were carried out at 2 different epochs. 

• Different reduction algorithms were used. 

• Different source detection thresholds were used. 

2.3.2 Region 2 

Two catalogues (Gregory et al. 1994; Wright et al. 1994) 
were produced with data from the PMN survey, Southern 
declination strip. Gregory et al.(1994) compare the results 
of two data reduction techniques used to compile the two 
source lists: 

• They find good agreement of positions of strong sources, 
but up to la disagreement in faint source positions. 

• They find the Wright et al.(1994) resulting catalogue 
contains more sources overall. 
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2.3.3 Problems within individual catalogues 

Two effects cause the flux fimit of the 87GB survey to vary 
with deciination, increasing from ~25mJy at 8 > 60° to 
~40mJy at 8 = 0° (Gregory & Condon 1991). One is the 
elevation dependence of the sensitivity of the 91m NRAO 
telescope, caused by a higher noise level at larger angles 
from the local zenith (8 = 38°). The other is the fact that 
adjacent observing tracks are spaced slightly further apart 
with decreasing declination. The Parkes telescope does not 
suffer decreased sensitivity at low elevations. However, the 
declination-dependent track spacing holds for low declina- 
tions (8 < —70°) and causes the flux limit to decrease from 
~50mJy near the zenith to ~40mJy near the pole. 

2.3.4 What we use as our catalogue 

We use the same cuts of the source catalogues as listed in 
Table 1 in Loan et al. (1997), resulting in sky coverage of 
~ 70%. We impose a working lower flux limit of 50 mjy. 
Figure |l| shows the distribution of the raw data over the 
sky, in equatorial coordinates. Here we discuss the procedure 
used to try correcting for the flux and number discrepancies 
mentioned in the discussion of the overlapping regions above. 
A measurement of A/" for each catalogue separately, at a 
given flux cutoff, shows that AfpMN-raw > Ns7GB-ra,w by 
a few percent for the range of flux cutoffs considered above 
as listed in columns 1-4 of Table [j]. It is possible that we 
are detecting actual large scale structure on the scale of 
the catalogues, but it is more likely that the inconsistencies 
(dicussed in Section |2|) in the methods by which catalogues 
were compiled are responsible. 

In an attempt to correct for the incompatibility of the 
two catalogues, we chose the following procedure to set 
Apa/jv ~ Ms7GB m - Mean density (A/ - ) matching - We take 
87GB fluxes as fiducial. For each flux cutoff (and corre- 
sponding M) in the 87GB ra m catalogue, we find the flux 
cutoff in the PMN ra m catalogue that gives Npmn as close 
as possible to N&tgb- Columns 1, 2, and 5-7 in Table [l] 
list these values. The numbers are not perfectly matched as 
the fluxes are quantizd in units of lmjy. The 60 mjy is the 
worst matched. 



3 SPHERICAL HARMONIC ANALYSIS 

It is common (e.g. Peebles 1980) to express the overdensity 
at a given point as a superposition of Fourier modes 8^: 

5(r) "7 (r)= (2^/ d3fc5ke " lkr ' (1) 

where r refers to comoving coordinates. In linear theory 
(<5(r) <C 1), the perturbation corresponding to a given 
wavenumber, k, is assumed to evolve independently of all 
other k modes. 

An alternative presentation of power on different scales 
is by Spherical Harmonic Analysis (SHA), which has become 
popular in recent years in studying the CMB and galaxy 
surveys. The idea of using SHA to study galaxy distribu- 
tions was first proposed over 20 years ago (Peebles 1973). 
The successful application of this method to actual data had 
to wait another decade for galaxy surveys (e.g. IRAS) that 



covered enough of the sky and contained enough objects to 
make SHA feasible. Unfortunately, for most of the objects 
detected, no redshifts were available. In this case, the true 
3D galaxy distribution is seen projected onto the surface of 
the celestial sphere and the spherical harmonic analysis is 
done in 2D (Scharf et al. 1992, Lahav 1994). When redshifts 
became available, the spherical harmonic analysis formalism 
was extended to 3D (Scharf & Lahav 1993; Fisher, Scharf 
& Lahav 1994; Fisher et al. 1995a; Heavens & Taylor 1995; 
Nusser & Davis 1994). In the 3D case, the redshift distortion 
introduced by individual galaxies' peculiar velocities must 
be corrected for in order to get distance information. As the 
radio catalogues described here provide only the angular po- 
sition of galaxies, we consider here 2D harmonic formalism, 
but relate it to the 3D power-spectrum. 

3.1 2D Harmonic Formalism 

Expanding a galaxy distribution in terms of spherical har- 
monic order, I, gives us information about clustering am- 
plitude as a function of angular scale, A9 ~ j, on the sky. 
This angular scale can be related to a linear scale by incorpo- 
rating distance information in the form of either individual 
galaxy redshifts or overall redshift distribution, N(z), of the 
galaxies. 

The distribution is expanded in spherical harmonics: 
a(f)=]Taryr(f), (2) 

/ . m 

where f is the angular position on the sky, and the resulting 
harmonic coefficients are 

JV 

aT = j2Yr\h). (3) 

In the case of incomplete sky coverage we derive the 
harmonic coefficients from the data as 

cz b s = Y ™* ( fi ) - N dn Y ™* ( 4 ) 

i=l Jn ots 

where M is the number of sources per steradian. The first 
term on the right hand side is aj™ from eqn. ^ and the sec- 
ond term is the correction for the unsurveyed and excluded 
regions of the sky. The second term (calculated by Monte 
Carlo integration) is included so that, for a Poisson distri- 
bution of galaxies, c™ obs = (for / > 0) when N ga i — > 00. 

3.2 Predicted Harmonics 

We can predict what the above af should be, for a given cos- 
mology and power spectrum of fluctuations, P(k) (dictated 
by a dark matter model). Note that the 2D spherical har- 
monic analysis formalism can be extended to other cosmolo- 
gies (Langlois et al, in preparation), but for this discussion 
we assume an Einstein-de Sitter universe (fi = 1, A = 0). 

We also account for the fact that a galaxy survey will 
only detect a fraction of the radio galaxies that populate the 
surveyed volume by including the selection function, <I?(r), 
which is the probability that a radio galaxy will be detected 
at distance r in this flux limited survey. The observed galaxy 
distribution is expanded in summation because galaxies are 
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c5 = + 9Q° 




<5 = -90° 



Figure 1. The distribution of radio sources (above flux cutoff of 70mjy in 87GB ra w, 74 mjy in PMN ram ) from the 87GB - PMN ra . 
catalogue areas as selected by Loan et al. 1997. This is an Aitoff projection in Equatorial coordinates. 



Table 1. Results of mean density (Af) matching procedure. The columns labeled % dif. give the percentage by which the larger Af 
exceeds the smaller Af. 



Slim Afg7GB-raw AfpMN-raw % dif. Sn m AfpMN-match % dif. 



50 


4457 


4929 


10.5 % 


54 


4433 


0.5 % 


60 


3520 


3836 


8.9 % 


61 


3744 


6.4 % 


70 


2870 


3123 


8.8 % 


71 


2891 


0.7 % 


SO 


2386 


2616 


9.6 % 


85 


2403 


0.7 % 


90 


2021 


2241 


10.8 % 


97 


2012 


0.4 % 


100 


1751 


1922 


9.8 % 


107 


1758 


0.4 % 


200 


655 


684 


4.4 % 


205 


658 


0.5 % 



discrete objects. To model the underlying mass distribution, 
which is continuous, the expansion is in integral form. Here 
we follow the formalism developed by Lahav, Piran & Treyer 
(1997) for the X-ray Background. As the sources of the X-ray 
Background are not resolved, the sources are flux weighted 
in their formalism. In our case, weighting by the flux of each 
radio source can make the measurements very noisy. Instead, 
we derive here the number weighted case (cf. Piran & Singh 
1996 for a similar formalism for Gamma-ray bursts). We 
begin by writing the harmonics as 

aT = J dV c {l + b R 6 m (r)]$(r) (5) 

and the fluctuations over the mean in the distribution (for 
I > 0) as 

aT = J dto dr r 2 S R (r) <E-(r) YT* (?), (6) 

where the comoving volume element dV c has been split into 
its radial and angular components, and <5r is the fluctuation 
in the number of radio galaxies with respect to the mean. 



For simplicity (reflecting our ignorance about the evolution 
of radio sources) we have assumed linear epoch-independent 
biasing factor between the radio sources and the mass fluc- 
tuations 

5 R (r) =b R S m (r). (7) 

Now we make use of the Fourier relation between the 
spatial and wavenumber density fluctuations. By substitut- 
ing the Rayleigh expression, 

e*' r = 4^^ji(kr)Yr*(f)yT(^), (8) 

I I'll 

into eqn. oL it becomes 




^(i ! )Tr(f) / d 3 ks k (z)Yr'(i)ji(kr). (9) 



Ira 

By substituting this expression for the density fluctua- 
tion term in eqn. [| we get 




dujdrr 2 $(r) ^ )* VT^Y?" (f) 



© 0000 RAS, MNRAS 000, 000-000 



Large scale structure in deep radio surveys 5 



x / d 3 k 8±{z)Y{T *(k) jl(kr). 



(10) 



By the orthogonality relation, J dujYr {^Y^" (r) = 5™, m ' 

(where 8™, m is the Kronecker delta function), the ^ over I 
and m drops out of the above expression. 

Mass density fluctuations evolve with time, growing as 
the scale factor 



5k(z) =<5k(0) (1 + z)- 



(11) 



In linear theory n = 1 for an Einstein-de Sitter universe, 
and from here on we adopt this value. 

We take the mean square of the a^'s and then the en- 
semble average of the result. We use the definition below to 
rewrite the density fluctuation term (<5k S^,) as: 



(«,) = (2^) 3 P(fc)«5 (3) (k-k'), 



(12) 



where P(fc) is the power spectrum of fluctuations in mass 
and 8^ is the 3-dimensional Dirac delta function. 

After grouping some of the terms into the window func- 
tion, *S?i(Je), we have 



<i«n 2 > = -bi / dfcfc 2 p(fc)i* ; (fc)! 2 



(13) 



The expression we evaluate now gives the ensemble av- 
erage (in the rms sense) harmonic coefficients. This is the 
prediction for the harmonic coefficients we would get if we 
averaged the a™'s measured by all possible observers (under 
the same observational conditions we have) in the universe. 
Of course, the radio data represent only one realization, but 
for each harmonic I, we have (21 + 1) independent measure- 
ments. 

To model the fluctuations we parametrize P(k) for a 
family of Cold Dark Matter (CDM) models. In this case the 
shape of P(k) is determined by the product of the density 
parameter and the Hubble constant Y = £lh. This parameter 
is related to the size of the horizon when the densities of 
matter and radiation where equal. For standard CDM F = 
0.5, while for low density CDM model which fits empirically 
clustering of APM and IRAS galaxies, V = 0.2. A smaller T 
gives more power on large scales. The normalization of the 
mass power spectrum is commonly given by as,M, the rms 
fluctuations in spheres of 8 h^ 1 Mpc. 

3.2.1 The Window Function, *;(fc) 

We have grouped the distance and redshift terms into "ifi(k) 



*i(fc) = / dr r 2 $>(r)ji(kr)(l + z\ 



(14) 



Here r is the comoving distance, which is related to redshift 
by r = — (1 + z)~ 1/2 ]. We do not know the selection 

function, &(r), explicitly and have to express in terms of 
the observed redshift distribution. Radio luminosity func- 
tions (Dunlop & Peacock 1990) can be used to calculate the 
number of radio sources we expect to detect as a function 
of redshift, z, given the frequency at which we are observing 
and a flux cutoff. The resulting redshift distribution, N(z), 
gives the number of radio galaxies observed per steradian at 
a redshift, z. 

N(z) was derived from the Dunlop & Peacock (1990) 
radio luminosity functions for a number of flux cutoffs, as 



explained in Foan et al.(1997). The luminosity range of radio 
galaxies is so large that the resulting redshift range probed 
is approximately the same for a range of flux cutoffs. 

From the definitions of 4>(r) and N(z), we know that 



dr r 2 $(r) =N= I dz BN(z) 



(15) 



where M is the number of galaxies per steradian in our sur- 
vey and B is a normalisation constant. Since the above equal- 
ity holds for all r, we can rewrite ^i(k) in terms of N(z). 



®i(k) = 2BQ / dx N(z(x))ji(x)(l - QxY 



(16) 



where we have set x — kr and Q = H /2kc. 

It is very insightful to plot the window function on top 
of the curves for fc 3 P(fc)/(27r 2 ) ~ (^) 2 , as it shows what 
scales of k are probed by particular statistic applied to a 
specific catalogue. The window function for the quadrupole 
|V>2(fc)| 2 (eqn. H for 70 mjy survey, is shown by the dashed- 
dotted line in Figure EL together with the power-spectra for 
standard CDM (with shape parameter T — 0.5) and low 
density CDM (r = 0.2) models. In both cases we assume 
normalization erg, m = 1 and &r = 1 (no biasing). The figure 
illustrates that radio surveys can serve as probes of scales 
fcjT 1 ~ 600/i _1 Mpc, between those accessible to galaxy sur- 
veys (e.g. APM) and the CMB (e.g. COBE). Other vari- 
ants of this plot which include other probes of the power- 
spectrum (e.g. local redshift surveys and the X-ray Back- 
ground) are given elsewhere (e.g. Baugh & Efstathiou 1993; 
Lahav et al. 1997; Wu, Fahav & Rees, in preparation). 



3.2.2 Shot Noise 

Finally, we need to consider the discreteness of the galaxies, 
which contributes a shot noise term to the data harmonics. 
To account for this, we add a shot noise estimate to the pre- 
dicted model harmonics. Each measured a™ will be subject 
to the same shot noise contamination 



W\ 2 ) tot 



i\ 771 I 2\ 

= (h I ) r 



. + 



\ a l I >SN- 



(17) 



(18) 



An added advantage of number weighted harmonics (over 
flux weighted) comes in here in the shot noise term. For the 
number weighted case the shot noise is finite, whereas for 
the flux weighted case it diverges if bright sources are not 
removed, (cf. Fahav et al. 1997). 



3.2.3 Incomplete Sky Coverage 

The radio surveys considered here only cover about 70% of 
the sky. The remaining 30% is filled in approximately uni- 
formly (according to eqn. |^) to the mean density of the sur- 
veyed sky. We do this to avoid the crosstalk between the har- 
monic coefficients, aj™, that can result from an incomplete 
sky. As we mention in Sections 5 and 6, we still measure 
significantly larger than predicted harmonics amplitudes in 
the data. However, the SHA formalism presented in Section 
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Figure 2. Radio surveys as probes of scales between those ac- 
cessible to galaxy surveys (e.g. APM) and the CMB (e.g. COBE). 
The solid and dashed lines represent k 3 P(k) / (2tt 2 ) ~ (^) 2 f° r 
standard CDM (r = 0.5) and low density CDM (r = 0.2) models, 
both normalized with <jg m = 1 at k ~ 0.15. The open squares 
at large fc's (small scales) are estimates of the power-spectrum 
from 3D inversion of the angular APM galaxy catalogue (Baugh 
& Efstathiou 1993, 1994). The elongated 'box' at small fc's rep- 
resents the COBE 4-yr CMB measurements (E. Gawiser, private 
communication). It corresponds to a quadrupole Q=18.0 fiK for 
a Harrison-Zeldovich mass power-spectrum, via the Sachs- Wolfe 
effect, or <Tg Af = 1.4 for a standard CDM model. The dashed- 
dotted line is the window- function \ip2{k)\ 2 for the quadrupole 
(eqn. [14] ) of radio sources brighter than 70 mjy and an Einstein- 
de Sitter universe. It illustrates the scales probes by the spherical 
harmonic analysis of radio surveys (eqn. hjj). It peaks at k^ 1 ~ 
6OO/1 -1 Mpc, indicating the capability of radio surveys to probe 
intermediate scales. The vertical scaling of the window function is 
arbitrary. The solid triangle and circle are estimates of the power- 
spectrum at fe*, assuming the shape of standard CDM and LD- 
CDM power-spectra respectively. They are based on the measured 
quadrupole from the combined 70 mjy 87GB — P MJV ma j c ), sam- 
ple, which correspond to crg n = buds^M ~ 9 and 5 for standard 
CDM and LDCDM respectively. Given the problems of catalogue 
matching and shot-noise, these points should be interpreted at 
best as 'representative upper limits', not as detections. 



3 is for full sky coverage. It is possible that despite our treat- 
ment of the unsurveyed regions, we are still seeing an effect 
of incomplete sky coverage. 

To make sure that we are comparing the data and model 
harmonics on equal footing, we want to deal with the unsur- 
veyed sky more rigorously. This can be done in two ways. 
One is to solve the inversion problem (given the known mask, 
i.e. unsurveyed region geometry) and to retrieve the data 
harmonics for a complete sky. The second approach is to 
adapt the model harmonics to additionally incorporate the 
effect of incomplete sky coverage and predict the rms har- 
monics for all observers with the same incomplete sky we 
have. We use the second method to examine the amplitude 
of the effect of incomplete sky coverage. 



Figure 3. Predicted ensemble average harmonics, using N(z) 
for a 70 mjy cutoff for LDCDM (higher curve) and CDM (lower 
curve) power spectra. Solid lines correspond to (l a ™| 2 |) mo£ j f° r 
full sky, and dotted lines correspond to (l c ; n | 2 |} mod f° r 70% sky 
coverage, for an ideal 4.85 GHz survey. 



The equation for the predicted harmonics for an incom- 
pletely surveyed sky is then 



l C i I )tot = z^ZJ^"' 



I )tof. 



(19) 



where the W^i m tensor models the unsurveyed regions (see 

Peebles 1973, Scharf et al 1992), and the (|a[/| 2 ) tot 's lire the 
predicted harmonics for a completely surveyed sky (equation 
[t|). This can also be written as {\c?\ 2 ) tat = (\cr\ 2 ) mod + 
(\ C T\ ) sn> wnere the masked shot-noise term is (|c™| ) SN = 
{n oba /4Tr)Af. 

Details of the calculation of the W^)" 1 tensor are given 
in Scharf et al (1992), and Zare (1987). The result of this 
correction is a lowering of the predicted ensemble average 
amplitudes, not an increase, as might have explained the 
large observed harmonic amplitudes. The corrected ampli- 
tudes (|Ci"| 2 ) mod in the absence of shot-noise for a 70 mjy 



J, compared to the predic- 
tion for full sky coverage (la™! 2 ) mod . This shows the trend 



flux cutoff, are shown in Figure | 

m I 2 \ 



in predicted harmonic amplitudes between the CDM and 
LDCDM power spectra and between full vs. partial sky cov- 
erage for a single flux cutoff. There is also a trend across 
chosen flux cutoff. With respect to the amplitudes shown 
for a 70 mjy cutoff, those for 50 mjy and 100 mjy are, 
repectively, ~ .5x and ~ 2x as high. 



4 THE PREDICTED DIPOLE 

The I = 1 harmonic term has additional significance in the 
study of large scale structure. The well known CMB dipole 
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anisotropy is thought to arise from our peculiar velocity rel- 
ative to the CMB. This motion is attributed to the gravita- 
tional pull of anisotropically distributed mass overdensities 
around us. In linear theory (Peebles 1980) the peculiar veloc- 
ity vector, v, is proportional to the gravitational acceleration 
vector, g, 

v oc — g, (20) 

where Q is the density parameter at the present epoch and 
b is the bias parameter (assuming linear biasing). 

Galaxy surveys do not represent the total mass directly, 
although the light emitted by galaxies at certain frequencies 
is likely to be indicative of their mass and thus related to 
the underlying mass distribution. The dipole can be calcu- 
lated as a flux weighted sum over each galaxy's position (see 
eqn. |33| for a more general form). If the galaxy's light is pro- 
portional to its mass this can represent the gravity vector as 
both the flux and gravitational force fall off with distance 
as 1/r 2 . 

In practice, however, all available galaxy catalogues are 
flux limited, and detect only a fraction of the galaxies in 
any volume surveyed. This means that the calculation will 
give an estimate of the dipole which is produced by only 
those galaxies detected out to a limited effective distance. 
Moreover, in this study, we have fluxes at radio wavelengths 
and the radio emission from galaxies may not be a good 
indicator of galaxy mass. We therefore prefer to calculate a 
number-weighted dipole. As a result, we cannot expect the 
number-weighted dipole to be aligned with the mass gravity 
vector and with the CMB Local Group velocity dipole. 

We now consider the various effects that may contribute 
to the dipole amplitude measured for a population of radio 
sources. We parameterise the dipole as 

^(0) = AT + A cos 6, (21) 

where the left side of the equation is the mean density of 
objects measured, per unit solid angle, as a function of di- 
rection, 6, on the sky. The angle 8 is measured from the 
direction of peak dipole effect. M — £L is the number of 
galaxies per steradian in the combined 87GB-PMN cata- 
logue for a given flux cutoff. A is the dipole amplitude we 
expect to measure due to the effects discussed below. The 
different notation for the dipoles is related as follows: 



A 



An 



-D 



\/<Kf>- 



(22) 



where D = | ^\ ?i| is the amplitude of the number- weighted 
dipole calculated and listed in Tables |§] and g, and (|ar| 2 ) 
is the dipole term of the harmonic prediction of our models 
(Section |). 

As described below, the velocity dipole is a separate ef- 
fect from the large scale structure dipole, although related 
to it. The large scale structure induces our motion, which in 
turn creates the blueshifted excess of galaxies in the direc- 
tion of motion. The measured dipole will be a combination 
of these two effects, plus the shot noise contribution. 

4.1 Velocity term 

For the sake of the following calculation, the population of 
radio galaxies is assumed to be isotropic in its rest frame. 



Ellis and Baldwin (1984) predicted the dipole amplitude we 
should measure in a flux limited radio galaxy survey, due 
to our motion with respect to the rest frame of the radio 
galaxy distribution. Consider a population of radio galaxies, 
whose radio continua can be characterized by the power law 
S oc v~ a ,a positive. In the direction of the sun's motion, 
galaxies are blueshifted, so the observed 4.85 GHz band ac- 
tually samples a lower frequency band in the rest-frame of 
the galaxy, and hence a higher flux. The overall effect is that 
we see more galaxies in the direction of motion and less in 
the opposite direction than we would if we were not moving 
with respect to the galaxies' rest frame. This deviation from 
isotropic number counts introduces a dipole pattern into the 
galaxy counts. This is similar to the 'Compton-Getting' ef- 
fect flux dipole due to the motion of the observer relative to 
a sea of radiation (e.g. the CMB or the X-ray Background). 

This velocity dipole effect was predicted by Ellis & 
Baldwin (1984) as 



Avel = [2 + x(l + a)]~ N, 



(23) 



where x is the slope of the integral source count per unit 
solid angle, above a given limiting flux, Su m , 



dQ {S > Slir ' 



oc S~ 



(24) 



As mentioned above, a is the slope of the power law spectra 
of the surveyed objects, and v is our velocity with respect 
to the frame in which the counts are isotropic. Note, that 
the above estimate is for a number weighted dipole. We take 
x = 1 ( = 1.5 in a Euclidean universe), and a = 0.75 (mean 
spectral index of radio galaxy spectra) ). The velocity in- 
ferred from the CMB dipole is v sun -CMB ~ 370^. The 
velocity dipole effect we should measure is 



Ar 



4.625 x 10" 3 A/\ 



(25) 



As an example of the magnitude of the effect, in the 
combined catalogue used below we have N = 2878 for 
Sum = 70mJy. In this case: 



A ve i — 13.3 



(26) 



In principle, we can subtract the velocity dipole based 
on the solar motion with respect to the CMB, and then con- 
sider the residual dipole as being purely due to large scale 
structure. It should be noted that the A ve i, unlike Alss and 
Asn discussed below, is not an rms quantitity with zero 
mean, but in fact a correctable value. It is thus quite dif- 
ferent in character from the other two terms. Eq. (23) only 
holds for a single population of radio sources characterized 
by spectral index a and counts index x. In reality, the 4.85 
Ghz sample is composed of different populations, and there- 
fore the expression for A vs i should be a superposition due to 
the different populations. At present the distribution func- 
tions for a and x at 4.85 GHz are poorly known, so it is 
difficult to evaluate the exact prediction for A ve i. 

At this point, it is important to emphasize the distinc- 
tion between the two reference frames with respect to which 
the CMB dipole amplitude and direction are defined. One 
is the heliocentric frame: v S un-CMB ~ 370fcm/s towards 
Galactic coordinates I — 264° ; b — 48° . As the measure- 
ments are done in the heliocentric frame, v SU n-CMB will 
determine the amplitude of A ve i. The other relevant refer- 
ence frame is the Local Group frame: vlg-cmb ~ 600fcra/s 



© 0000 RAS, MNRAS 000, 000-000 



8 A. Baleisis, 0. Lahav, A. J. Loan, and J. V. Wall 



toward I = 268° ; b = 27° . Since the radio galaxies in the 
catalogues are outside the Local Group, this pertains to the 
LSS dipole prediction in our discussion, but as explained 
above there is no simple relation between a number-weighted 
dipole and the CMB Local Group velocity dipole. 

4.2 Large Scale Structure term 

This is the / = 1 term, from the predicted spherical harmon- 
ics in Section [| 

As an example of the magnitude of the effect we have 
assumed a power spectrum which fits the observed local 
clustering of optical and IRAS galaxies, parameterised as 
an unbiased, low density CDM (hereafter LDCDM) power- 
spectrum with shape parameter F = 0.2. We also assume an 
Einstein-de Sitter universe (Q = 1) and hence that fluctu- 
ations grow as (1 + z) -1 - We adopt a redshift distribution, 
N(z), corresponding to a flux cutoff of 70m Jy. For these pa- 
rameters the rms prediction, eqn. |D| gives (|a™| 2 } = 83.52, 
which translates into 



Alss = 7.7. 



4.3 Shot Noise term 



(27) 



As mentioned earlier, a shot noise term comes in because we 
are using discrete galaxies as tracers of a continuous quantity 
(fluctuations in the mass density) . We calculated the Poisson 
shot noise contribution to be 
1 ■ 

47T ■ 



(\ a T\ 2 )sN 



1 = AT, 



(28) 



A S N = 



Af. 



(29) 



In the case of our combined catalogue for a 70mJy flux 
cutoff, Af = 2878 and we expect: 



45.4. 



(30) 



Unfortunately, the shot noise turns out be the dominant 
component of any dipole we measure in the 87GB-PMN sur- 
vey. It is simple to estimate the number of objects we would 
need to detect the velocity and large scale structure dipoles 
at the same level as shot noise. For the velocity dipole, we 
equate eqns. M and f29 and find that we need 



Af ~ 3.35 x 10 4 galaxies per str, 



(31) 



which translates into N ~ 4 x 10 J galaxies over the sky. The 
large scale structure dipole, Alss is proportional to Af (see 
eqns. [l3| [l5|, |l6|), so by equating eqn. |2^ and the expression 
for the large scale structure dipole we find we need 

Af » 1.06 x 10 5 galaxies per str, (32) 

(or TV ~ 1.3 x 10 6 galaxies total) to be able to detect a large 
scale dipole (for the LDCDM power spectrum, the shape of 
N(z) and the other parameters mentioned above). 

The predicted velocity, large scale structure, and shot- 
noise dipole effects (for the example parameters mentioned 
above) are summarized in Table g. Our a priori calculation 
shows that we are unlikely to detect a dipole (and higher 
harmonics, see below) in the 87GB-PMN catalogues due to 
the high level of shot-noise expected. 



Table 2. The amplitudes of various dipole effects, compared 
to the number weighted dipole we detect (after the Af matching 
procedure) in the 87GB — PMN ma t c h catalogues, for a 70mjy 
flux cutoff. 



Dipole 


A (amplitude) 


Velocity 


13.3 


Large Scale Structure 


7.7 


Shot Noise 


45.4 


Detected 


73.3 



Nevertheless, we have attempted to calculate the dipole 
and higher harmonics from the catalogues for several rea- 
sons: (i) The detection of the angular correlation function 
(Loan et al. 1997) and the Supergalactic Plane (Section 7) in 
87GB raw and PMN raw suggests the existence of large scale 
structure, (ii) It may well be that our formalism of Gaus- 
sian random fields does not fully characterize real structure, 
(iii) Comparison of the predictions to observations can serve 
as a consistency check of the validity of the 87GB ra m and 
PMN r aw flux calibration, (iv) The application to existing 
catalogues can serve as a pilot study for future deeper sur- 
veys. 

As described below, we find that the observed dipole is 
much higher than expected. We attribute this discrepancy to 
the flux matching of the catalogues. However, it is interesting 
that after applying an ad-hoc flux matching procedure we 
still detect a dipole (see Table ^ larger than predicted. Flux 
cutoff's of 50mJy and lOOmJy result in qualitatively similar 
results - the detected dipole is still higher than the predicted 
shot noise, large scale structure and velocity dipoles. The 
magnitude of this surplus of the detected dipole with respect 
to the predicted dipoles increases with decreasing flux limit. 
We discuss the implications for this detection in the next 
few sections. 



5 DIPOLE CALCULATION USING THE 
COMBINED 87GB - PMN match 
CATALOGUE 

Despite the numerous uncertainties in relating the radia- 
tion emitted by galaxies to the underlying mass distribu- 
tion, angular dipole calculations for a variety of galaxy cat- 
alogues (e.g. Meiksin & Davis 1986; Yahil, Walker & Rowan- 
Robinson 1986; Villumsen & Strauss 1987; Lahav 1987; 
Harmon, Lahav & Meurs 1987; Lahav, Rowan- Robinson 

6 Lynden-Bell 1988; Plionis 1988; Lynden-Bell, Lahav & 
Burstein 1989; Kaiser & Lahav 1989; Scharf et al. 1992), 
have found amazingly good alignment (10° — 30°) to the 
CMB dipole in the Local Group frame. Similar alignment 
has been detected using IRAS redshift surveys (Rowan- 
Robinson et al. 1991; Strauss et al. 1992; Webster, Lahav 
& Fisher 1997) and optical redshift surveys (e.g. Hudson 
1993). 

The detection of a dipole closely aligned with the CMB 
dipole, in a number of surveys at various wavelengths, sug- 
gests that it is worth calculating the dipole for the radio 
galaxies in the 87GB-PMN survey. We only have angular 
coordinates and fluxes. The 87GB-PMN surveys have a me- 
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Table 3. Initial dipole calculation for 87GB — PM N raw . Smi n 
is the flux cutoff in mjy, M is the number of galaxies per steradian, 
D is the amplitude of the dipole, and I and b give the Galactic 
coordinate direction in which the dipole is pointing. 



Cumulative dipole (vii = 1) 





M 


D 


1 


b 


50 


4671 


1342.3 


299° 


-21° 


70 


2985 


709.2 


297° 


-15° 


100 


1831 


468.4 


300° 


-17° 



dian redshift of z ~ 1. This is important for the question 
of convergence of the measured galaxy dipole to the CMB 
dipole. Above some large enough distance from us, we ex- 
pect the distribution of galaxies to be homogeneous, and not 
contribute to the CMB dipole. The 87GB-PMN catalogues 
are deep enough to probe the convergence of the dipole far 
beyond the limit set by IRAS and optical surveys. 

On the other hand, we do not know how well radio radi- 
ation from a galaxy indicates its mass, or how radio galaxies 
are biased with respect to the underlying mass distribution. 
In addition, the sparse sampling of the 87GB-PMN surveys 
may make it difficult to even detect a dipole signal. If we 
assume an Einstein-de Sitter cosmology, for a 70mJy cutoff 
there are only ~ 500 galaxies out to 100 h^ 1 Mpc and ~ 1500 
objects out to z ~ 0.1 in the combined radio catalogues. 

After filling in the unsurveyed regions of the sky, we 
calculated the radio galaxy dipole for a range of flux cutoffs 
in the survey. The general form of the dipole calculation (for 
a full sky distribution of N galaxies) is 

JV 

D = ^Wifi. (33) 

i=l 

For a flux weighted dipole, Wi = fi, the flux of galaxy i. 
We use Wi = 1 , and calculate the number weighted dipole 
because, although the flux weighted dipole is the optimal 
calculation to recover the acceleration vector, it is not al- 
ways reliable. As mentioned above, only some wavelengths 
of a galaxy's radiation are good indicators of its mass. If 
there is no correlation between L ga i at a given wavelength 
and Mg a i , the flux weighted dipole will not give an estimate 
of the gravitational pull we are experiencing. In such a case, 
it is better to use a strictly number weighted dipole. Note 
that this no longer gives an estimate of g because it does 
not incorporate the 1/r 2 fall off in flux to each galaxy (as 
with Wi — fi), and as such, is no longer directly comparable 
to the LG-CMB dipole. The results of our initial calculation 
are shown in Table |j. As mentioned, no alignment was ex- 
pected. The measured amplitudes, however, are significantly 
larger than the number weighted prediction from Section 3. 

Note that the above dipole calculations have been done 
under the assumption that the 87GB raw and PMN raw cat- 
alogues were produced using similar reduction criteria, and 
that for each flux cutoff Npmn- raw = NsiGB-raw (where 
M is the number of galaxies per steradian on the sky). As 
mentioned in Section [], this is not a valid assumption. For 
the combination of the 87GB raw and PMN raw catalogues 
that we use, we have found the regions of the sky covered by 



Table 4. Recalculated dipole - S m in is the flux cutoff in 
mjy, Af is the average number of galaxies per steradian in the 
87GB — PMN matc i x catalogues, D is the amplitude of the dipole, 
and I and b give the direction in which the dipole is pointing in 
Galactic coordinate. 



Cumulative dipole (wi = 1) 





M 


D 


1 


b 


50 


4443 


415.3 


284° 


25° 


60 


3622 


742.7 


293° 


-17° 


70 


2878 


307.2 


286° 


15° 


80 


2392 


251.7 


282° 


16° 


90 


2016 


169.2 


287° 


28° 


100 


1753 


158.4 


298° 


16° 


200 


656 


45.3 


274° 


-19° 



the PMN r aw catalogue to have larger M values for a given 
flux cutoff. 

The resulting dipole in the flux matched catalogue, for a 
70mJy flux cutoff was D — 307.2, a significant decrease from 
the original value. The direction also changed considerably, 
to I = 286°, b = 15°, 21° away from the LG-CMB dipole, 
38° away from the Sun-CMB dipole. 

It is interesting to note that our derived dipole (Table 
is closer to to the Local Group CMB dipole than to the 
heliocentric CMB dipole. As previously mentioned, however, 
not much can be inferred from the alignment between our 
number-weighted dipole and the LG-CMB dipole. At most, 
this suggests that the velocity dipole does not dominate the 
measured dipole. 

The decrease in the measured dipole amplitude was en- 
couraging enough to recalculate the dipole for other flux 
cutoffs. The results for a range of flux cutoffs are shown in 
Table ^ and Figures ^| and ^. The dipole amplitudes are 
smaller by a factor of 2 — 3, and show much closer align- 
ment to the CMB dipole than before. For the 60 mjy cutoff, 
the dipole amplitude is still quite high. We attribute this 
to the poor flux matching for this flux cutoff which left a 
number density mismatch comparable to some of the origi- 
nal A/pMiv-ra™, A/g7Gs-raiD percentage differences at other 
flux cutoffs. 

We see a significant reduction in the amplitude and 
change in the direction of the dipole, after applying fairly 
naive (and seemingly small) changes to the catalogues. It is 
possible that the catalogue geometry (which we have tried 
to correct for by approximately uniform filling in of the un- 
surveyed regions) is causing a residual dipole signal. 

5.1 Summary and discussion of dipole calculation 

For the radio sources in the 4.85 GHz, 87GB-PMN combined 
catalogue, for a flux cutoff of 70 mjy, we find the following: 

(i) We predict the effects of our motion and of large scale 
structure to be comparable in amplitude to within a factor 
of 2. 

(ii) We expect shot noise to dominate any dipole in the 
87GB-PMN catalogues. 

(iii) When we measure the dipole we find that, 

• For no correction of the catalogues, A is 3.5cr above 
the noise, but over 50° away from the LG-CMB dipole. 
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Figure 4. After A/" matching : Detected dipole amplitude, D, for 
flux cutoffs from 50 mjy to 200 mjy. 



■~- IRAS l.gjy 



260 280 300 

Galactic longitude, 1 



Figure 5. After Af matching : Detected dipole direction, for 
flux cutoffs from 50 mjy to 200 mjy (filled circles), the IRAS 1.2 
Jy dipole (open pentagon), LG-CMB dipole (open square), and 
Sun-CMB dipole (open triangle). 



• When we try the flux matching procedure, we mea- 
sure a dipole signal of only 0.5cr above noise, but now 
the dipole points 20° away from the LG-CMB dipole. 

We attribute most of the initially measured dipole sig- 
nal to the fact that two different reduction algorithms were 
used in compiling the 87GB-PMN catalogues. The signal 
decreases after flux matching, but is still above shot noise 
level. In addition, we now see closer alignment with the LG- 



CMB dipole, that we should not detect a priori. What might 
produce these effects? 

We hope that, after the flux matching done, the 
87GB — PMN match catalogue is consistent, and not pro- 
ducing additional significant dipole effects, resulting from 
the slight, remaining N mismatch. However, as mentioned, 
there is the possibility that survey geometry is significantly 
influencing the dipole alignment. 

An alternate explanation of the dipole detection over 
what we predict, could come from the large scale structure 
prediction. The large scale structure dipole we predict is 
for an ensemble average of observers over a universe (with 
LDCDM power spectrum). The dipole measurement of the 
87GB-PMN catalogue represents only one realization, which 
may greatly differ from the ensemble average. We may live in 
a highly non-homogeneous part of the universe. There is the 
possibility that we are detecting real large scale structure. 
As we report in Section 7 we detect a slight increase in 
the mean density of sources in the 87GB raw and PMN raw 
catalogues within 10° — 15° along the Super Galactic Plane. 
This structure may contribute to the remaining dipole signal 
we measure. 

Some support for a dipole pattern in the distribution of 
radio galaxies comes from an analysis by Lahav & Shaver 
(1991, unpublished). They analysed an optically-radio se- 
lected sample (Shaver 1991) of 92 galaxies which flux at 408 
MHz larger than 1 Jy and blue magnitude brighter than 
14.5, 2/3 of them with redshift z < 0.02. If the dipole is 
so local, it is due to only few hundreds of sources (see our 
discussion of the SGP below). 



6 SHA APPLICATION 

We now compare observations and predictions for higher 
harmonics. The fluctuations in radio sources can be related 
to the fluctuations in mass via the linear bias parameter ba, 
as as, r = bi?o"8,M- The models we use for our predictions are 
for a universe with f2 = 1 geometry, but we allow 2 possi- 
ble power-spectra: (i) standard CDM power spectrum with 
Q — 1.0, h — 0.5, and normalization specified by the rms 
fluctuations in radio counts in 8 h~ l Mpc spheres, (ts,r = 1.0 
and (ii) a power spectrum parameterised as an unbiased 
ip%,R = 1.0) low density CDM (LDCDM) spectrum shape 
parameter V = 0.2 (Bardeen et al. 1986). 

We can rewrite the corrected harmonic coefficients as 



cr 



A/-(fi oi)s /47r)(2Z + l)^ 



|2 

\tot 



1 + 



<icrr>. 



■(34) 



for |c™| tot being the predicted harmonics (signal+noise, for 
an incomplete sky) or the observed harmonics. Subsequently, 
Ci tells us directly how many standard deviations above the 
shot noise a measurement or prediction is. The predicted 
harmonic in the absence of noise (|c™| 2 ) mo( ; and Ci (includ- 
ing noise) for a flux cutoff of 70mJy and both power spectra 
are listed in Table |E| Incomplete sky coverage has been 
corrected for, according to eqn. 

The models predict (see Table |B|) that due to our ob- 
serving constraints, we will not see much power (if large 
scale structure exists in CDM or LDCDM form) over shot 
noise for any order I. The shot noise level is at Ci = 1. 
However, we have already measured a larger than predicted 
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Table 5. The predicted harmonic coefficients for CDM and LD- 
CDM models for a flux limit of 70 mjy, for I = 1 to I = 5. 





CDM 




LDCDM 


1 


1 ) mod 


Ci 


(l^l 1 )mod 


Ci 


1 


16.4 


1.0040 


50.7 


1.0120 


2 


17.3 


1.0042 


47.3 


1.0116 


3 


17.5 


1.0043 


43.6 


1.0106 


4 


17.0 


1.0042 


39.6 


1.0096 


5 


16.5 


1.0040 


36.4 


1.0088 



dipole amplitude in the 87GB — PMN match data, so we are 
interested in what the data harmonics will show. 

Before N matching, for the 87GB-PMN catalogue, we 
see a large dipole (I = 1) and significant signal in harmonics 
up to I = 10. The 87GB - PMN ma tch data show ampli- 
tudes somewhat reduced (mainly for the lower harmonics 
/ < 6, but they are still far above the predicted amplitudes. 
Even for a full sky, filled to the same mean density as for the 
70mJy flux matched case, the harmonic coefficients oscillate 
considerably around the shot noise prediction. It is also in- 
teresting to note that the flux matched data harmonics have 
higher overall amplitudes the lower the flux cutoff, following 
the general trend predicted in the models. 

Figure ^ shows 

the harmonics for the 87GB — PMN match catalogue with a 
70mJy flux cutoff, compared to the predicted harmonics for 
the same flux cutoff, corrected for the incomplete sky cov- 
erage. For the model prediction, we plot the resulting C;'s 
for (\ c r\ 2 )mod x 1, Xl5, x30, and x45, and can see that the 
galaxy harmonics are consistent with being about a factor 
of 15 to 30 above the prediction (including shot noise). 

The above predictions for the radio harmonics assume 
no biasing, i.e. the normalization is fixed to U8,r = ft/jo^M = 
1. If we interpret the observed, N matched harmonics as 
being real, we can use them to estimate the amplitude o"8,k 
for a given shape of the power-spectrum. For example, for 
the quadrupole (I — 2) and flux limit of 70 mjy the ob- 
servation C2 ~ 1.3 can be made to agree with the CDM 
and LDCDM predictions if <J8 : r ~ 9 and ~ 5 respectively. 
(With better data one should estimate the normalisation us- 
ing Maximum Likelihood, e.g. Scharf et al. 1992; Fisher et al. 
1994.) These crude results suggest a very high bias parame- 
ter 6_r, but in qualitative agreement with the measurement 
of high correlation length for radio sources. For example, 
for 87GB and PMN at redshift 2 = 0, f(r) = (r/r )" 18 
with ro = I8/1 -1 Mpc (Loan et al. 1997), assuming sta- 
ble clustering. This can be translated (Peebles 1980) to 

cr S ,R « >/ L86 ( r o/8) 1 - 8 ~ 2.8. 

Figure ^ shows (solid triangle and circle) the estimates 
of the power-spectrum at fc _1 ~ 6OO/1 -1 Mpc, assuming the 
shape of standard CDM and LDCDM power-spectra respec- 
tively. They are based on the measured quadrupole from 
the combined 70 mjy 87GB — PMN match sample. Given 
the problems of catalogue matching and shot-noise, these 
points should be interpreted at best as 'representative up- 
per limits', not as detections. 
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Figure 6. C; comparison for 87GB — PMN matc h 70mjy data 
and LDCDM model predictions {{\cY l \ 2 ) mod x 1, Xl5, x30 and 
x45) The dots connected by a solid line indicate the data har- 
monics, while the various dashed and dotted lines indicate the 
model harmonics. C; = 1 is the shot noise level. 

7 THE SUPERGALACTIC PLANE (SGP) 

The large-area uniformity of the 87GB-PMN surveys was 
adopted as the criterion for a sensible choice for sky mask 
and flux-density limit. This produced the masked catalogue 
used to calculate the dipole and higher harmonics in this 
study, as well as the angular correlation function in Loan 
et al. (1997). Nevertheless, nearby structures, such as the 
Supergalactic Plane (SGP), may introduce detectable inho- 
mogeneities into the catalogues on large angular scales. 

The SGP is an over-dense region of the local universe. 
William Herschel was the first to note that 'nebulae' seemed 
to be more concentrated in a band across the heavens. More 
recently, the Virgo cluster was seen as the centre of a 'meta- 
galaxy' or 'Local Supercluster'. De Vaucouleurs (1975) rec- 
ognized a great circle along which there was an increased 
density of galaxies in the Shapley-Ames optical galaxy cat- 
alogue. This great circle defines a system of Supergalac- 
tic spherical coordinates, in which the SGP lies along the 
equator and the north pole lies at Galactic coordinates 
(£ = 47.37°, b = 6.32°). The SGP is readily apparent as 
a region of greatly enhanced surface density of galaxies in 
more recent large-area surveys of nearby optical galaxies ( 
e.g. Lynden-Bell & Lahav 1988). Although the SGP seems 
to be a local phenomenon, perhaps analogous to the sheets 
and filaments discovered by galaxy surveys, it is not clear 
where it ends. Studies of the SGP using the IRAS 1.2 Jy 
and ORS surveys (Lahav et al., in preparation) show that 
the SGP is not a simple planar structure. 

7.1 Histograms approach 

In the case of full sky coverage, the great circle along which 
the density of galaxies is enhanced can be found by calculat- 
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Figure 7. Source surface density of 87GB raw, as a function of 
SGB and flux-density limit. Error bars are Poissonian. The 9 
bins are evenly spaced in sin(SGB). The area of each bin is es- 
timated by Monte Carlo integration. The dotted line shows the 
mean source density over the whole survey. 



Figure 8. Source surface density of PMNraw, as a function of 
SGB and flux-density limit. Error bars are Poissonian. The 9 
bins are evenly spaced in sin(SGB). The area of each bin is es- 
timated by Monte Carlo integration. The dotted line shows the 
mean source density over the whole survey. 



ing the 'moment of inertia' (Webster, Lahav & Fisher 1997; 
Lahav et al., in preparation). This approach is obviously 
only valid for a uniformly-selected catalogue with full sky 
coverage, which unfortunately is not the case here. More- 
over, clusters and voids outside the 'plane' may confuse the 
analysis as well. We prefer instead a more direct approach, 
where each catalogue (87GB ra w and PMN raw ) is analysed 
independently by histograms of number counts. 

In order to explore whether the SGP is indeed vis- 
ible in the 87GB ra w and PMN raw catalogues, we now 
calculate the surface density of sources as a function of 
position in Equatorial (a, 8), Galactic (£,b) and Super- 
galactic (SGL, SGB) coordinate systems. The &7GB raw and 
PMNraw catalogues are treated here separately, to avoid 
being mislead by the flux-density miscalibration highlighted 
earlier. In the study of the SGP we use the more conservative 
flux limit of 100 mjy. 

In Supergalactic coordinates, the sky is divided into 
strips of roughly equal area, equally spaced in sin (SGB) and 
covering the full range of SGL. The density of sources in each 
strip is found by adding up the number of sources in each 
catalogue within each strip, and dividing by the masked area 
of the strip (that is, the area of the strip inside the mask 
for each catalogue). The total area of a strip may be easily 
calculated using spherical geometry, but this simple calcu- 
lation is complicated by the numerous parts of the sky that 
are excluded from the analysis by the choice of the mask. In- 
stead, the area of each strip within the mask is estimated by 
Monte Carlo integration. Many points are put down at ran- 
dom positions within each strip, and the number of points 
included within the masked area is calculated. Thus, the area 
of each strip is the ratio of included random points to the 
total number of random points, multiplied by the total (un- 
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Figure 9. Surface density of sources above 100 mjy in the 
87GB raw and PMN ra w catalogues, as a function of SGB. Error 
bars are Poissonian. The 15 bins are evenly spaced in sin(S'GB). 
The area of each bin is estimated by Monte Carlo integration. Re- 
sults for each catalogue are calculated separately, and normalized 
by the average number density in that catalogue. 
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masked) area of the strip. This estimate of the area is found 
to be accurate to within a small fraction of 1 per cent. 

Figures ^ and ^ show how the number of sources in the 
87GB raw and PMN raw catalogues changes with the flux- 
density limit and as a function of SGB. Both catalogues 
appear to be complete above 50mJy, but display incom- 
pleteness at the 35mJy limit (not shown). The 87GB ram 
plot, Figure [?], is consistent with a uniform distribution at 
all flux limits, although there is a suggestion of a small ex- 
cess near SGB = 0° . The PMN raw plot, Figure |, displays 
a strong density enhancement near the Supergalactic equa- 
tor (SGB — 0°) at all complete flux limits. 

Figure |^ shows how the surface density of sources above 
a flux-density limit of 100 mjy varies as a function of SGB 
in the 87GB ra w and PMN ra m catalogues. There is a sig- 
nificant increase in source surface density close to the SGP 
in the PMN raw catalogue, but any increase is not statis- 
tically significant in the 87GB raw catalogue. Reploting this 
histogram for a different binning scheme retained this ex- 
cess around the SGP, indicating that these result are not 
sensitive to the choice of bin size. Tables ^| and ^ quan- 
tify the excess density within the central strip (either side 
of SGP = 0°) for the 87GB ra w and PMN raw catalogues, 
respectively. Approximately 70% of the SGP in each hemi- 
sphere is included in the central strip. Repeating this anal- 
ysis with strips of finite width in SGL, both for the whole 
range of SGB and also for sources within 10° of the SGP, 
shows that no single cluster of sources causes this increase, 
but that the over-density close to the SGP is distributed 
along the Supergalactic equator. 

A similar procedure was followed in Equatorial coordi- 
nates, plotting similar histograms with bins equally spaced 
in sin(<5) and covering the whole range of a It is clear that 
87GB raw and PMN r aw have slightly different flux calibra- 
tions, in accord with the discussion earlier. Even far above 
the published flux-density limits of both surveys (~ 50 mjy) 
there are more sources per unit area above a given flux limit 
in PMN raw than there are in 87GB raw- Fortunately, since 
both surveys are analysed separately, this problem will not 
affect these results. However, it does mean that the flux 
limits given here are somewhat different from the real flux 
limits. This difference will not greatly affect the generally 
observed trend, because the radio luminosity function is so 
broad. There is also the possibility of biases caused by ig- 
noring sources which lie just below the flux limit through 
experimental error in the surveys. Apart from the disturb- 
ing flux mismatch highlighted here, no significant density 
enhancements are seen in Equatorial coordinates. 

No significant patches of increased surface density are 
seen when a similar analysis is performed in Galactic coor- 
dinates (£,b), plotting histograms of bins equally spaced in 
sin(fe). Failure to detect any increase in source surface den- 
sity close to the Galactic Plane shows that there is very little 
contamination from Galactic sources. 

Similar analyses were conducted, binning the sources in 
bins of equal width in right ascension («), Galactic longitude 
(£) and Supergalactic longitude (SGB). Again, no significant 
departures from the mean sources density were observed. 



7.2 Comparison with previous studies 

The PMN r aw survey shows a significant enhancement of 
the radio source surface density within ~ 10° of the Su- 
pergalactic Plane: an increase of ~ 15% compared to the 
mean value. This is good evidence for concentration of ra- 
dio sources towards the SGP, although it is puzzling that 
there is no corresponding over-density in the northern sky. 
The concentration of low-Supergalactic-latitude sources in 
PMNraw is not due to one cluster of sources, but rather 
it is spread along the whole SGP. The number of sources 
necessary to create this enhancement is quite small ( ~200). 

The tests for increased numbers of radio sources close to 
the SGP follows a similar analysis of the Molonglo 408 MHz 
survey of the southern hemisphere (Shaver & Pierre 1989). 
Shaver & Pierre found a significant increase in the surface 
density of sources close to the SGP. This work confirms the 
presence of significant concentration of radio sources within 
10° of the standard SGP equator in the southern hemi- 
sphere (PMNraw), but this is not mirrored in the north 
(87GB ram ). Shaver (1991) pointed out that there may be an 
asymmetry in the SGL-distribution of the radio galaxies (in 
the Great Attractor direction), and this would explain the 
fact that the SG concentration is more evident in the south- 
ern hemisphere than in the north. Shaver & Pierre (1989) 
also analysed a small redshift survey, finding that there was 
a significant preference for nearby (z < 0.02) radio sources 
to lie close to the SGP, but that this preference was not ob- 
served for radio sources at higher redshift (0.02 < z < 0.1). 

The present analysis finds an over-density of about 130 
radio sources with a 4.85 GHz flux density above 100 mjy 
within an area of about 1500 square degrees, up to 10° either 
side of the SGP in the southern hemisphere. A less signifi- 
cant over-density is also observed in the north. The models 
of Dunlop & Peacock (1990) predict that the surface density 
of radio sources above 100 mjy within a redshift z < 0.02 is 
~18 per str (1 str=3282. 97 square degrees). It is interesting 
to note that the over-density found along the SGP in the 
PMNraw survey could be accounted for by approximately 
the same number of radio sources found within a redshift 
z < 0.02 (~ 60h~ 1 Mpc). 

If we assume that the over-density in the PMN raw 
survey is caused by nearby radio galaxies out to a redshift 
z = 0.02, then we can estimate the thickness of SGP. The an- 
gle subtended by the over-dense strip is 6 ~ 20°, and thus, 
at the distance corresponding to z = 0.02 (d — cz/Ho = 
60/i~ 1 Mpc), the SGP has a thickness ~ 2dtan(<9/2) « 
20h~ Mpc. This also raises the question of how far the SGP 
extends. One test of the hypothesis that the over-density 
close to the SGP really is caused by nearby radio sources 
would be to cross-correlate the radio sources with those in 
the optical IRAS catalogues, in a similar way to the study 
by Shaver (1991). 

It may be possible to cross-identify objects that are 
both optical IRAS sources and also radio sources in the 
87GB T aw and PMN ra w catalogues. Repeating this analy- 
sis after removing radio sources that appear in the optical 
and IRAS catalogues may show whether the apparent over- 
density close to the SGP is caused by a flattened distribution 
of nearby sources. 

It could be argued that these results may be pure sta- 
tistical fluke, and that the detected over-density is caused 
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Table 6. Excess number of radio sources close to the SGP 
for 87GB raw , above 100 mjy. Total number = 8272, Area 
= 15732 deg 2 , mean density = 0.526 deg -2 . The number ex- 
pected (N) is calculated from the mean density. The distance 
from the mean is the number of standard deviations (\/N) 
that the detected source excess (N — N) lies from the expected 
number. Note that the mean density is an over-estimate since 
it includes any over-density close to the SGP. Thus, the dis- 
tance from the mean is under-estimated. 



Central bin Bin area # found # expected # excess Distance from mean 



width (deg) 


(deg 2 ) 


(AT) 


(AT) 


(N - N) 


(N 


±3.82 


1103 


602 


580 


22 


0.9 


±6.38 


1819 


989 


956 


33 


1.0 


±11.54 


3237 


1741 


1702 


39 


0.9 


±19.47 


5335 


2828 


2805 


23 


0.4 



- n/Vn) 



Table 7. Excess number of radio sources close to the SGP for 
PMN r aw, above 100 mjy. Total number = 7829, total area = 
13486 deg 2 , mean density = 0.581 dcg~ 2 . 



Central bin 
width (deg) 


Bin area 
(deg 2 ) 


# found 

(AO 


# expected 
(A 7 ) 


# excess 
(N - N) 


Distance from mean 
(N - N/y/N) 


±3.82 


914 


613 


531 


82 


3.3 


±6.38 


1534 


1020 


891 


129 


4.1 


±11.54 


2766 


1742 


1606 


136 


3.3 


±19.47 


4604 


2807 


2673 


134 


2.5 



by a relatively small number of radio sources. However, it is 
not clear how a bias towards detecting sources that lie close 
to the SGP could have entered the radio surveys. 

On large angular scales, there is a clear detection of 
the SGP in the southern sky, but not in the north. The 
over-density in the south may be attributed to nearby radio 
galaxies lying preferentially toward the SGP, but it is not 
clear why there is no significant detection in the north. Pos- 
sible explanations of this conundrum include the possibilities 
that nearby sources are somehow excluded from 87GB raw , 
or that they are counted more than once in PMN ra w ( e.g. 
nearby double- lobed radio sources may be included twice). 
Cross-identification of sources in the radio surveys with op- 
tical and IRAS galaxies (both relatively shallow surveys, 
out to m 6000 km/s) could be used to find the nearby ra- 
dio sources, to check if they are indeed responsible for the 
over-density close to the SGP. 



8 DISCUSSION 

We have studied in this paper the possibility of using radio 
surveys with median redshift z ~ 1 to probe large scale 
structure. In particular, we have utilised the technique of 
spherical harmonics to predict the amplitude of the dipole 
and higher in the angular distribution of radio galaxies. 
Our conclusions are: 

(a) The dipole is due to 2 effects which turn out to be 
of comparable magnitude: (i) our motion with respect to the 
CMB, and (ii) large scale structure, 

(b) Catalogues like 87GB and PMN have the potential 



of probing structure on large scales, e.g. The quadrupole 
measures scales of ~ 600fe _1 Mpc. 

(c) Unfortunately, the Poisson shot noise in these sparse 
catalogues is expected to be large than the clustering signals, 
as predicted for a family of Cold Dark Matter models. 

(d) However, we detect dipole and higher harmonics 
in the combined 87GB — PMN raw catalogue which are far 
larger than expected. We attribute this to a 2% flux mis- 
match between the two catalogues. Ad-hoc corrections to 
match the catalogues may suggest a marginal detection of a 
dipole. 

There are two likely explanations for this. One is that 
we are still seeing the effect of survey geometry, or incompat- 
ibility of data reduction algorithms in producing the cata- 
logues. The second is that, because the catalogues represent 
only one realization of the ensemble average harmonics we 
predict, we may be detecting local structure that causes the 
large amplitude we measure. 

(e) To detect a dipole and higher harmonics unambigu- 
ously, a catalogue with full sky coverage and ~ 10 6 sources 
is required. 

(f) We have investigated the existence and extent of 
the Supergalactic Plane in the above catalogues. In a strip 
of ±10° of the standard Supergalactic equator we find a 3cr 
detection in PMN raw, but only la in 87GB raw . This analy- 
sis demonstrates that large-scale structure studies based on 
new radio surveys nearing completion will provide important 
estimates of distant structure and its evolution. On-going 
surveys with the VLA, FIRST (Becker, White & Helfand 
1995) and NVSS (Condon 1997), will yield approximately 
10 6 sources over the sky, thus overcoming the limitations 
imposed by the shot-noise at the current 4.85 GHz surveys. 
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Furthermore, the positional accuracy of these new surveys 
will by sufficient for optical spectroscopy using candidates 
directly from the catalogues. Although difficult to obtain, 
redshifts to the radio sources are essential in order to im- 
prove the measurement of the power-spectrum by localizing 
the clustering pattern in 3 dimensions. 
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